METHOD AND APPARATUS USING PSEUDO-INVERSES OF 
LINEAR TRANSFORMATIONS IN MULTI-CARRIER 
MODULATION RECEIVERS AND TRANSCEIVERS 

Technical field 

5 This invention relates to the use of pseudo-inverses and truncated pseudo- 
inverses of linear transformations to create time progressions received by 
Multi-Carrier Modulation (MCM) receivers, in particular in Orthogonal 
Frequency Division Multiplex (OFDM) receivers, including but not limited to, 
IEEE 802.11 compliant receivers and transceivers, as well as their use to 

10 determine availability of service for one or more MCM communications 
protocols. 

Background Art 

Figure 1 depicts a prior art receiver of a Multi-Carrier Modulated (MCM) 
signal. The MCM signal is sensed by antenna element 100 and fed 102 to 
is module 110. Module 110 contains a Low Noise Amplifier (LNA) 114 and 
possibly some form of switch diagrammatically depicted as 116 which 
collectively feed 112 down converter 120. 

The transmitted frequency may use any of several frequency bands, 900 
MHz, 2-2.5 GHz and 5-6 GHz, being commonly used transmission bands. 

20 Down converter 120 generates a down converted version of the sensed 
signal, which is fed 122 to Variable Gain Amplifier 124, which generates an 
amplified, converter signal which is fed 126 to Band Pass Filter 130. Band 
Pass Filter 130 removes undesirable noise components introduced by the 
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down conversion and amplification, generating a filtered intermediate 
frequency signal 132. 

The filtered intermediate frequency signal 132 is coherently split into two 
coherent intermediate frequency signals 134 and 136. Intermediate 

5 frequency signal 134 is presented to mixer 160 and intermediate frequency 
signal 136 is presented to mixer 170. Local Oscillator (LO) 140 generates a 
stable reference signal 142 which is split into two coherent reference signals 
144 and 146. Reference signal 144 is presented to mixer 160. Reference 
signal 146 is presented to phase offset 150, which generates a phase offset 

]0 re f erence signal 152 which is presented to mixer 170. Phase offset 150 
imparts the equivalent of a phase shift of 90 degrees onto phase offset 
reference signal 152 with respect to reference signal 146. 

Mixer 160 uses intermediate frequency signal 134 and reference signal 144 to 
create a first intermediate frequency component signal 162 in a frequency 
15 range compatible with A/D converter 190. 

Mixer 170 uses intermediate frequency signal 136 and phase offset reference 
signal 152 to create a second intermediate frequency component signal 172 
in a frequency range compatible with A/D converter 210. 

Signals 162 and 172 may include undesirable frequency components 
20 requiring further filtering before presentation to A/D converters 190 and 210, 
respectively, but these filters have not been shown in the interests of clarity. 

Often the first and second intermediate frequency component signals 162 and 
172 contain signals in a frequency range under 100 MHz, in many cases on 
the order of 1-20 MHz. The maximum frequency range is often known as the 
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intermediate frequency and determines the band pass frequency range of 
band pass filter 130. 

A/D converters 190 and 210 respectively generate first sample data stream 
192 and second sample data stream 212. The sampling rate is usually a 
5 multiple of the intermediate frequency, which by the Nyquist theorem is 
theoretically 2, and often in practice at least 2.5. 

The sampled data streams 192 and 212, also known as the I and Q sample 
streams, are sent to digital processor 250. In many situations, they are 
merged, buffered 260, then conveyed 262 across a communication network 
io 270, and delivered 272 to processing engine 280. 

The sample sizes of the A/D converters vary for specific applications, but may 
be any of 6, 7, 8, 9, 10, 12 or more bits per sample. These digitized samples 
are often packed into a computer format of 8, 16, or 32 bits. 

Processing engine 280 will often place these merged digitized samples into at 
15 least one sample input buffer 300 residing in memory 290. Memory 290 is 
accessibly coupled 282 to processing engine 280. Based upon these 
buffered, digitized samples 300, processing engine 280 will often utilize a 
heuristic mechanism to determine whether the sampled channel sensed by 
antenna 100, and processed by the above discussed mechanism is active, 
20 generating a Clear Channel Assessment (CCA) flag 310. 

The prior art has focused upon Clear Channel Assessment based upon 
classical signal detection theory. Such developments focus Bayesian 
detection mechanisms, which calculate the probability of successful detection 
against the probability of false detection. Such mechanisms tend to require 
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long start-up times, preferably capturing at least one start of burst or message 
header containing a training sequence, allowing timing synchronization 
between the user and the transmitting communications medium. 

A/D converters 190 and 210 usually sample their respective input signals at 
5 the same rate, in fact most often sampling those signals using a carefully 
constructed clocking scheme controlling clocking skew between them. 

One or more clocks may control the timing of processing engine 280. One 
clock may be operating above 100 MHz, and perhaps operating at much 
higher rates, such as 240 MHz or higher. 

10 Processing engine 280 may transfer the merged data samples from 
temporary buffer 260 based upon the triggering of an interrupt, in some cases 
using a Direct Memory Access (DMA) mechanism (not shown). The DMA 
mechanism may operate across network 270 to transfer the digitized input 
samples to sample inputs 300 residing in memory 290. Note that memory 

15 290 may include both volatile and non-volatile memory components. 

One or more program counters may control the operations of processing 
engine 280. Processing engine 280 fetches one instruction for each of the 
program counters to control the operations of processing engine 280. Note 
that a single program counter would be compatible with a Single Instruction 
20 processing engine, whereas multiple program counters would be compatible 
with a Multiple Instruction processing engine. 

Alternatively, the operation of processing engine 280 may be controlled by a 
collection of states, such as the one-hot state machines often found in FPGA- 
based designs. As used herein, a program step will refer to at least an 
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instruction or processor control state providing the controls necessary to 
execute one or more steps of the inventive method. A program system as 
used herein will refer to the collection of program steps implementing an 
embodiment of the inventive method. 

5 Inverse linear transform 340, as found in the prior art, approximates the 
inverse of a non-singular linear transform which was used to create a signal 
progression. That signal progression transports across the physical transport 
layer(s) of the communication protocol to create the transported version of the 
signal progression received as the sample list in buffer 260 and subsequently 

10 found in sample inputs 300. 

As used herein, MCM refers to the communication of a data stream by 
dividing that data stream into multiple parallel sub-streams, each having a 
lower bit rate, and then concurrently modulating these sub-streams with 
separate sub-carriers. The separate carriers may or may not be isolated from 
is each other. Frequency Division Multiple Access (FDMA) protocols, including 
AMPS and GSM, use frequency bins separated by guard bands as the sub- 
carriers as in Figure 2A. Such MCM protocols require steep bandpass filters 
that completely separate the sub-carriers. 

Orthogonal Frequency Division Multiplexing (OFDM) uses densely spaced 
20 sub-carriers and overlapping spectra as shown in Figure 2B, eliminating the 
need for steep bandpass filters. OFDM sub-carriers are, by construction, 
mutually orthogonal within the protocol specified sampling window. In many 
cases, both the transmitter and receiver employ complementary Fast Fourier 
Transforms (FFT) and Inverse Fast Fourier Transforms (1FFT) to transmit and 
25 receive the data sub-streams. 
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OFDM has been studied for use, in conjunction with Direct Sequence Spread 
Spectrum techniques, to create CDMA-OFDM protocols as shown in Figure 
2C. The transmitter would first apply a Waish-Hadamard transform to the 
multiple data sub-streams to create multiple spread data sub-streams. The 
5 multiple spread data sub-streams would then be transformed by an IFFT to 
create an intermediate frequency modulated signal, which is then up- 
converted to the transmission frequency band. 

The receiver would down convert 120 (and filter 130) the amplified antenna 
reception 112 to create a received intermediate frequency 132 as discussed 
10 in Figure 1 above. This received intermediate frequency signal 132 would be 
split (143 and 136) and mixed (160 and 170, respectively) with a local 
reference 144 and a phase offset version of the reference 152, which would 
be sampled by A/D converters 190 and 210, respectively. The output of A/D 
converters 190 and 210 are the I and Q sample streams as discussed above. 

is OFDM protocol research has lead to the specification and deployment of 
communications protocols in a variety of application areas including, but not 
limited to, Digital Video Broadcast, Digital Audio Broadcast, and wireless data 
networks. 

The following formulae provide a first of two equivalent definitions of Walsh- 
20 Hadamard transforms as used in spread spectrum communications. 
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The H matrices use an alphabet of two symbols, 0 and 1. In such an 
alphabet, the complement of 0 is 1 and the complement of 1 is 0. The over 
bar marks in formulae (2) to (4) refer to taking the component-wise 
5 complement of each element of the matrix involved under the bar. 

Formula (1) depicts the matrix, which is a 1 by 1 matrix. Formula (2) 
depicts H 2 , the 2 by 2 matrix generated as shown from H v Formula (3) depicts 
H 4 , the 4 by 4 matrix generated as shown from H 2 . Formula (4) depicts H 2N the 
2N by 2N matrix generated from H N , where N is a power of two. 

10 The following formulae provide a second equivalent definition of Walsh- 
Hadamard transforms as used in spread spectrum communications. 
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Note that various developments of the G matrices may involve a normalization 
factor, which has not been included. The G matrices use an equivalent 
alphabet of two symbols, -1 and 1, for which the complement of -1 is 1 and 
5 the complement of 1 is -1 . The over bar marks in formulae (6) to (8) refer to 
taking the component-wise complement of each element of the matrix 
involved under the bar. 

Note that the G matrices are more often used in practice, because the 
absolute value of every entry in the G matrices is the same. The IS-95 
10 communications protocol defines 64 logical channels encoded by G 64 . 

The H matrices are often used for pedagogical purposes or as part of a 
process leading to the G matrices, providing a more accessible relationship 
with the standard definitions of bits. 
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As used herein a bit represents an alphabet possessing two symbols. 
Multiple bits are the concatenation of single bits, preferably representing a 
single alphabet. 

The prior art also includes Discrete Wavelet Transform (DWT) coding, which 
5 is a powerful extension of the linear transform coding discussed to this point. 

The following formulae define Discrete Wavelet Transforms as found in the 
prior art. 



ft 0 ff° 

" - i f 'o 


4 


4 
















ft ,n~l 


of 1 


DO 








mt" S) 











1 



m 



(10) 



at , % *a% , t mb 3 * 



(II) 



Where the Kronecker deltas are defined as 

if h - s\ then 6*-*' - 1 elbe - 0 
if I — I s ^ Www S IM = 1 else hi i? ~ 0 

Formula (9) shows a matrix A with m rows and an unlimited number of 
columns. A is defined as a wavelet matrix of rank m if formulae (10) and (11) 
are satisfied by the components of A, a s k . These components usually belong 
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to an algebraic sub-field of the complex numbers, such as rational complex 
numbers, real numbers, rational real numbers, or the complex numbers 
themselves. To simplify the discussion, the components of A will be assumed 
to be complex numbers. The over bar of formula (11) refers to taking the 
complex conjugate of the expression under that bar. 
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Formula (12) defines a sub-block matrix A, of matrix A containing the columns 
of A from l*m to l*m+m-1. Formula (13) shows a second way of looking at A 
as an arbitrary long vector with entries A,. Suppose that only finitely many of 
the A, components are non-zero. Further suppose that A N1 is the first non- 
zero component and A N2 is the last non-zero component. Let g=N2-N1 +1 . 



Formulae (14) and (15) define the Laurent series A(z) for the matrix A defined 
by formula (9) in two different ways. 
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From hereon, this discussion will assume that only finitely many of the A, 
components are non-zero. Such matrices A will be known as discrete wavelet 
matrices. To further simplify the discussion, from hereon N1 will be taken to 
5 be zero. As one of skill in the art will realize that N1 could be non-zero and 
the resulting matrices would essentially be equivalent to "shifted" matrices 
where N1=0. 

Any DWT is equivalent to a matrix A, which has m rows and m*g columns, 
where m is called the rank and g is the genus of the transform. Each DWT 
io transform is further characterized by an m*m characteristic Haar matrix which 
equals A(1). When g is one, the transform has a square matrix equal to its 
Haar characteristic matrix. 

Formula (11) asserts that the rows of a wavelet matrix a s =(a s 0 , a 5 ^^) have 
length m 1/2 and that they are pair-wise orthogonal when shifted by an arbitrary 
is multiple of m. 

The vector a 0 is often called the scaling vector, low pass filter, or scaling filter. 
The vectors a s for 0<s<m are often called wavelet vectors, high-pass filters, or 
wavelet filters. Formula (10) states that the sum of the components of the 



scaling vector is m, whereas for each wavelet vector, the sum of components 
is zero. 



The rank m of a transform corresponds to the sampling rate and to the 
number of bands in an m-band filter bank implementation. If a filter has rank 
5 m, then it samples the signal m times per unit time. When m is infinite, the 
sampling is continuous and the filter is analog. 

The genus g of the transform represents the number of symbols or signaling 
intervals over which the filter operates. When the rank m is finite, m*g equals 
the number of taps in each sub band filter. Note that if g is infinite, the filter 
10 has infinite duration and is not practicable. 

Discussion of the characteristic Haar matrix would entail a digression, which is 
not central to the invention. Suffice it to say that for a given rank m, the 
choice of Haar characteristic matrices ranges over a continuous (m-1) A 2 
dimensional family of matrices. 

15 Increasing the rank of m corresponds to increasing the spectral resolution 
whereas increasing the genus of a transform corresponds to increasing the 
overlap of successive transform windows. 

Note that Fourier transforms as well as Walsh-Hadamard transforms both can 
be defined, extended and analyzed by DWT techniques, though this is a topic 
20 well outside the range of this invention. 

Walsh-Hadamard matrices are a special case of Hadamard matrices. 
Hadamard matrices are square matrices of rank m containing components 
whose values are either +1 or -1. Hadamard matrices further satisfy 
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H T H=HH T =ml. Note that a Hadamard matrix may be an n*n matrix of rank m, 
where m is less than n. 

Consider a linear transform from a domain N dimensional space to a range, 
which is a second N dimensional space. If the linear transform does not cover 

5 the range, in other words, if there is a point in the second N dimensional 
space for which there is no point in the domain which transforms to that point, 
then the linear transform is singular. If for every point of the second N 
dimensional space, there is exactly one corresponding point in the first N 
dimensional space which transforms to that point, then the transform is non- 

10 singular. 

Linear transforms are well known to have matrices associated with them. 
When the linear transform is non-singular, its matrix is non-singular, has non- 
zero determinant and possesses an inverse, which is also non-singular with 
non-zero determinant. 

15 When the linear transform is singular, its matrix is singular, has determinant 0 
and does not possess a non-singular inverse. 

A further problem occurs when a linear transform goes from an N dimensional 
space to an M dimensional space, where N and M are distinct. Again, there is 
a matrix associated with the transform, but the transform cannot be non- 
20 singular. 

Work by a number of people, including Moore, Penrose, and Drazin, has led 
to a theory of inverses applicable to singular matrices, giving rise to several 
sometimes distinct, pseudo-inverses of a matrix. As may be expected, the 



pseudo-inverse of a non-singular (necessarily square) matrix is the classic 
inverse, which is non-singular. 

As used herein, R(A) will refer to the range of the linear transform for the 
associated matrix A and N(A) will refer to the null space of the linear 

5 transform for the associated matrix A. The addition of two vector spaces over 
the same field (which will usually be C, the complex numbers) is the vector 
space including exactly all linear combinations of the vectors of the two vector 
spaces. Note that much of this discussion is applicable to vector spaces over 
algebraic fields in general and in specific, almost always to algebraic sub- 

10 fields of the complex numbers including the rational real numbers, real 
numbers, and rational complex numbers, as well as the complex numbers. 
The discussion from hereon will focus on vector spaces over the complex 
number field for convenience and is not intended to restrict the scope of the 
claims herein. 

15 The conjugate transpose of a matrix A will be denoted herein as A* and will 
include components for a given row and column which are complex 
conjugates of the component of the column and row of A. 

Given a matrix A of m rows and n columns of rn*n complex number 
components, a matrix G of n rows and m columns of n*m complex number 
20 components is called herein an (i,j,k)-inverse of A if G satisfies the ith, jth and 
kth Penrose conditions: 

1. AGA = A 

2. GAG-G 
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3. (AG)* = AG 

4. (GA)* = GA 

The set of all (l,j,k)-inverses for A will be denoted by A{l,j,k}. 

The following are some basic facts about some of the various classes of 
5 inverses developed in detail in Generalized Inverses of Linear 
Transformations by S. L. Campbell and C. D. Meyer, Jr., © 1979, first 
published by Dover in 1991, ISBN 0-486-66693-X, particularly in chapter 6. 

G belongs to A{1} if and only if Qb is a solution of Ax^b, for every vector b in 
the range of A. This type of inverse is denoted (1)-inverse and is called the 
io Equation Solving Inverse. 

If G belongs to A{1,2}, then N(A)+R(G)=C n and R(A)+N(G)=C m . Each (1,2)- 
inverse defines complementary subspaces for N(A) as well as R(A). 
Conversely, for each pair of subspaces (P,Q), where P and Q are 
complementary to N(A) and R(A), respectively, uniquely determine a (1,2)- 
15 inverse, G PQ with R(G PQ )=P and N(G P Q )=Q. This type of inverse will be called 
a prescribed range/null space inverse. 

Two vector subspaces of a vector space will be referred to as complementary 
if the only element they have in common is the origin, and if linear 
combinations of vectors from these two subspaces cover the whole vector 
20 space. 

G belongs to A{1,3} if and only if Gb is a least squares solution of Ax=b for 
every vector b in C m . Gb will be referred to as a least squares solution of Ax=b 



when the distance between the hyperplane Ax and the vector b is minimal at 
Gb. This type of inverse will be called a least squares inverse. 

G belongs to A{1 ,4} if and only if Gb is the minimum norm solution of Ax=b for 
every vector b in R(A). This type of inverse will be called a minimum norm 
5 inverse. 

As used herein, the norm of a vector b is formed as the square root of the 
product b and b*. Note that the product of b and b* is a non-negative real 
number. The minimum norm solution Gb has the least norm of any solution of 
Ax=b. 

10 A{1,2,3,4} contains exactly one element, denoted as A + herein. A + is the 
(R(A*), N(A*))-inverse for A. A + b is the minimal norm least squares solution of 
Ax=b for any b in C m . If b belongs to R(A), then A + b is the minimal norm 
solution of Ax=b. A + is known elsewhere as the Moore-Penrose Inverse. 

Computing the Moore-Penrose inverse A + from the above definitions involves 
15 an unpleasant fact. If A is of neither full row rank nor full column rank, then the 
rank of A may be perturbed in an arbitrarily small way, dramatically changing 
the value of A + (see page 247 of Campbell and Meyer for a discussion and 
proof). 

Another approach to calculating the Moore-Penrose inverse A + involves use of 
20 the Singular Value Decomposition Theorem (see pages 6, and 247-262 of 
Campbell and Meyer). Matrix A is factored into A=UEV, where U and V are 

unitary (square) matrices and E has the form 

E = [ Diag(Eigen(A*A) 1/2 ) 0 ] 
[0 0 ] 
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The E matrix includes the positive eigenvalues of the matrix of the square 
roots of A*A. 

Defining A + = V*E + U*, can be shown to provide a well defined pseudo-inverse 
varying continuously with A regarding suitably chosen matrix norms (pages 
247-262 of Campbell and Meyer, matrix norms are defined and discussed on 
pages 210-213). This definition is often taken as the standard for these 
reasons. 

While much more can be said about this topic, the above definition is 
computationally demanding. Matrix inverses of non-singular square matrices 
using Gaussian elimination take on the order of N 3 operations. Calculation of 
all the eigenvalues of a matrix (A*A) 1/2 is an even greater task. A number of 
specialized algorithms are discussed in Campbell and Meyer, as well as in 
other literature sources which are much faster and often useful, but lack the 
generality of the Singular Value Decomposition derived A + . 

Note that the Moore-Penrose inverse, as well as {i j,k}-inverses, provide either 
some form of solution, or least-squares solution, for a linear algebraic system. 
While possessing many important qualities, these inverses lack some other 
desirable qualities. Let A and B be n*n complex matrices, there is no class 
C(i,j,k) of {i,j,k}-inverses A' and B~ for A and B respectively which imply any of 
the following: 

1. A A" = A" A, 

2. (A~) p = (A p )~ for all positive integers p, 

3. X is an eigenvalue of A if and only if 1M, is an eigenvalue of A\ 
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4. A P+1 A~ = A p , for positive integers p. 

The Drazin inverse and group inverse have at least these properties. 

Before discussing the Drazin inverse, we need to define the index of linear 
transformation A from C n to C n , which will be denoted herein as lnd(A). Ind(A) 
5 is the smallest non-negative integer k such that C n = R(A k ) + N(A k ), where A 2 
= A applied to A, A m+1 = A applied to A m , for any positive integer m. Note that 
if A is non-singular, Ind(A) = 0 and that Ind(O) = 1 . Further, if k = Ind (A), then 
R(A k )=R(A k+1 ). 

There are two ways to define the Drazin inverse of a square matrix A 
10 associated with linear transformation A from C n to C n . 

Let A be a linear transformation on C n such that k = Ind(A). Let A n =A 
restricted to R(A k ). Let x=u+v belong to C n , where u belongs to R(A k ) and v 
belongs to N(A k ). is invertible and define A D x = A^ 1 u. A D is the Drazin 
inverse of linear transformation A. This definition is known as the Functional 
15 Definition of the Drazin inverse. 

The Algebraic Definition of the Drazin inverse defines A D in C n * n for A in C n * n 
with lnd(A)=k as a matrix satisfying the following: 

A D A A D = A D 

A A D = A° A and 

20 A k+1 A D = A k . 

These definitions are equivalent, and for any A in C n * n , A D exists and is unique. 
If A is non-singular, then A D is exactly the standard matrix inverse. 
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Note that A° is not always a {1}-inverse for A, it doesn't always solve Ax=b. In 
fact A D b is a solution of Ax=b if and only if b belongs to R(A k ), where k = 
Ind(A). 

In one important case, when Ind(A)<=1 , A D solves Ax=b. In such cases A D is 
5 called the Group Inverse of A and is often denoted by A # . When it exists, A # 
can be alternatively defined as the unique matrix satisfying 

A A # A = A 

A*A A # = A # 
AA # = A # A 

10 Finally, the relationship between the Moore-Penrose inverse A + and the 
Drazin inverse A D for a matrix A in C n * n , A + = A D if and only if A + A=A A + . 

Modern radio receivers such as depicted in Figure 1 often face situations 
where the active reception channels to be decoded are a subset of all the 
channels supported by the protocol. Note that these channels are the 
15 received version of signals transmitted after being linear transformed by any 
of a number of matrices. These transform matrices may or may not be square 
matrices. They may not be invertible, even if they are square matrices. 

An example of this is found in the IEEE 802.11a protocol, where there are 
always 12 null frequency bins, but during the header transmission there are 
20 only 12 active data frequency bins out of the 64 frequency bins in the protocol. 
While only a small part of the frequency bins are required, there is no 
technique available to describe and control DSP resources at any level finer 
than a linear transformation. As a consequence, vital computational resources 
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are expended when only a small part of those resources need be used. 
During the header, only 12 out of the 64 results of the FFT are used. Note that 
the range of the encoding transformation for the preamble has a dimension of 
12 out of 64. The range of the encoding transformation of the body has a 
5 dimension of 52 out of 64. 

Figure 2D depicts an example emitted power spectrum requirement for 
transmitted OFDM signals as found in Figure 120 of the IEEE 802.11a 
protocol specification. 

"The transmitted spectrum shall have a 0 dBr (dB relative to the maximum 
10 spectral density of the signal) bandwidth not exceeding 1 8 MHz, -20 dBr at 1 1 
MHz frequency offset, -28 dBr at 20 MHz frequency offset and -40 dBr at 30 
MHz frequency offset and above. The transmitted spectral density of the 
transmitted signal shall fall within the spectral mask, as shown in Figure 120. 
The measurements shall be made using a 100 kHz resolution bandwidth and 
3 5 a 30 kHz video bandwidth." (17.3.9.2 Transmit spectrum mask page 28 of the 
IEEE Standard 802.1 1a-1 999 document) 

What is needed is a method of specifying the use of DSP resources based 
upon inverses of these matrices, which may not necessarily be square 
matrices, and which may be singular matrices even when they are square 
20 matrices. What is further needed is a way to determine the type of inverse of 
these matrices, as well as calculate the matrix inverse type, most useful for 
the specific radio reception problem. 

Additionally, it is desirable for software radios like that shown in Figure 1 be 
able to receive communications in multiple communications protocols. These 
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communications protocols may use the same frequency range for dramatically 
different protocols. An example of this is the use of the AMPs frequency 
channels by the IS-95 communications protocol. 

1S-95 employs a physical transport layer made of either a single or pair of 
physical channels reusing the AMPs physical channels. An IS-95 physical 
channel takes up 41 contiguous AMPs physical channels. When a single IS- 
95 physical channel is implemented, it is surrounded by a guard band of 9 
AMPs physical channels on either side in the frequency domain. When dual 
IS-95 channels are implemented, the pair of IS-95 physical channels are 
immediately adjacent to each other, with 9 AMPs physical channels on either 
side of the pair of IS-95 physical channels. 

What is needed is a mechanism to specify linear transformations suitable for 
rapidly reconfiguring reception by the software radio receiver between distinct 
communications such as AMPs and IS-95, allowing the reception of channels 
based upon suitable matrix inverses of these diverse linear transformations. 

Cellular telephone users regularly complain about coverage limitations in the 
United States. A given area will often support some of the wireless protocol 
standards, which may include AMPs, GSM and IS-95, but not all of them. 
Most cellular telephones today are built around transceivers, which 
communicate using only one standard. Cellular telephone transceivers need 
to sense which protocols are actively supported in the area near which the 
transceiver is operating. 

Cellular telephone users who travel internationally face a very similar problem. 
Again, most cellular telephones respond to only one of the common standards 
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of today, and there are several distinct standards employed in large areas of 
the world. These cellular telephones are often useless in areas not 
supporting that one standard for which they are compatible. Cellular 
telephone transceivers again need to sense which protocols are actively 
5 supported in the area near which the transceiver is operating, 

Summary of the invention 

Aspects of the invention address at least each of the above-mentioned needs. 

Many communications protocols involve a collection of communication 
channels collectively forming the dimensions of a finite dimensional vector 
10 space, of which, at any point in time, only a subset of those channels or 
dimensions must be received. The encoding of messages onto these 
channels is a time progression in the actively used dimensions of the vector 
space. 

Examples of this include the OFDM protocol IEEE 802.1 1a, where there are 
is always at least 12 of the 64 frequency bins unused in generating the time 
domain sequence which is transmitted. Another example is the GSM cellular 
telephone protocol, which uses a limited number of time slots and frequency 
bins to communicate necessary systems control and timing information to any 
transceiver in its service area for the initiation of a phone call in the service 
20 area of a base station. Similarly, IS-95 reserves certain logical channels, 
which are dimensions in the Walsh-Hadamard matrix-derived vector space for 
similar communication of necessary systems control and timing information to 
user transceivers attempting phone call initiation in the service area of a base 
station. 
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The invention utilizes a truncated version of the pseudo-inverse for the 
actively used portion of the transform to process received sample lists into at 
least one signal parameter. The actively used portion of the transform may 
not be all that the external transmitter is sending, but preferably is the actively 
5 used portion that the receiver wishes to receive. 

The inverse of the actively used portion of the transform is a pseudo-inverse, 
and the invention may employ distinct forms of pseudo-inverses even for the 
same communications protocol. The truncated version of the pseudo-inverse 
may generate a result including at least some of the active channels. The 
10 truncated version of the pseudo-inverse may also generate exactly the active, 
subspace results. 

The invention may use a truncated version of the actively used portion of the 
transform augmented by at least one signal which will be cancelled in the 
receiver to process received sample lists into at least one signal parameter. 

is By way of example, let a chirp signal s be added to a data-bearing OFDM 
time sequence symbol involving the 52 actively used subcarriers. On the 
transmitter side, the intermediate frequency output will be the product of a 
matrix M by the sequence of time stimulus vectors X of 53 dimensions, 52 
data bearing dimensions and one chirp signal dimension. The matrix MM will 

20 contain the 64 by 52 DFT matrix M concatenated with a 64-eiement chirp 
vector s. MM=[M s] has the array dimensions of 64 by 53. Note that the chirp 
signal s is required to comply with any communications protocol requirements 
such as emitted power spectrum constraints. 
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On the receiver side, compute N as the pseudo-inverse of MM, having array 
dimensions of 53 rows by 64 columns. Let NT be the matrix formed by 
truncating the row associated with s in N, which in this example is the last 
row. Estimate X = NT * [Received Samples]. 

Certain preferred embodiments of the invention support the reception of 
OFDM signals compliant with the IEEE 802.11a standard by processing a 
received sample list using a truncated pseudo-inverse of the trimmed matrix 
associated with the 64 point FFT (and/or DFT). It is further preferable to use 
at least two truncated pseudo-inverses of the trimmed square matrix, i.e. 
rectangular matrix associated with a trimmed FFT and/or DFT, the first 
excludes 12 null channel-frequency bins and the second excludes 52 null 
channel-frequency bins. The first is used during the data portion of a 
transmission and the second is used during reception of a transmission 
header. 

Note that by knowing what signal is being cancelled, NT can be pre-computed 
and stored. In the case of data-bearing IEEE 802.11a OFDM signals, up to 
64-52 = 12 signals may be actively cancelled. 

Additionally, when an IEEE 802.11a receiver wakes up after quiescence, it is 
preferable to use a third truncated pseudo-inverse of the FFT generating a 
minimal number of frequency bin samples. These frequency bin samples are 
then used to calculate a first energy estimate and a second energy estimate 
for the purposes of determining the CCA_flag. The CCA_fiag is set to busy 
whenever the first energy estimate exceeds the second energy estimate 
multiplied by a threshold. Note that at least some of the minimal frequency 
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bin samples may in effect integrate more than one of the frequency bins 
specified by the standard. 

Certain preferred embodiments of the invention include receivers, therefore 
transceivers, capable of supporting at least two communications protocols in 
wireless communications applications including cellular telephones and 
personal digital assistants. When such receivers wake up from quiescence, 
they process at least one sample list received from an electromagnetic 
transponder using a truncated pseudo-inverse allowing them to determine 
whether there is service support for at least one of the supported protocols of 
the receiver. 

It is preferable that if two communications protocols share the same frequency 
window, as in the case of AMPs and IS-95, that the truncated pseudo-inverse 
contains a third truncated pseudo-inverse and a fourth truncated pseudo- 
inverse. Preferably, the result of applying third and fourth truncated pseudo- 
inverses aids in determining whether there is service support for the first 
wireless protocol and second wireless protocol sharing the frequency window, 
respectively. This minimizes latency in determining which communications 
protocols are locally supported. 

If two communications protocols supported by the receiver/transceiver have 
separated enough frequency bands requiring separate down-conversions, 
then it is preferable in certain embodiments of the invention that separate 
sample lists be received for each protocol. These separate sample lists may 
be further processed using at least partially differing truncated pseudo- 
inverses to determine the availability of service for the different protocols. 
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The receiver is coupled to at least one electromagnetic receptor, receiving a 
sample list of at least two digitized samples based upon the electromagnetic 
receptor. Often these digitized samples are generated as the output of at 
least one A/D converter. The electromagnetic receptor may include at least 
5 one antenna element, and receiving the sample list may be derived from an 
electromagnetic field proximate to the antenna element. The electromagnetic 
receptor may include at least one semiconductor element, and receiving the 
sample list may be derived from an electromagnetic field based upon the bulk 
transport properties of the semiconductor receptor element. 

io The receiver may be coupled to multiple electromagnetic receptors. In such 
cases, receiving the sample list may preferably include receiving the sample 
list from at least two of the electromagnetic sensors. Examples of such 
embodiments include transceivers coupled to multiple directional antennas. 

Similarly, base station receivers must find each transmitting user's signal 
15 against the perpetual reception of many forms of background noise. Base 
stations may know where in a first vector space of reception to look for a 
user's signal, which may include but is not limited to, frequency slots, TDMA 
time-frequency slots, and spread spectrum time-aligned spreading 
sequences. 

20 There is a second vector space, a base station or other receiver must contend 
with, regarding the relative location of a transmitter. Often such receivers are 
coupled to multiple antenna elements and/or electromagnetic receptors, 
providing sensor streams from multiple overlapping reception lobes. The 
propagation effects from a transmitter to and/or from these antenna elements 

25 and/or electromagnetic receptors can be modeled as a linear transform of the 
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time progression of encoded signals at the transmitter derived from the lobe 
plot of the radio network. These linear transforms are almost never singular. 

What is needed is a mechanism for processing the received sample lists from 
multiple electromagnetic receptors and/or multiple antenna elements within at 
5 least one electromagnetic receptor to improve either location resolution for 
receiving a transmitted signal and/or improve the signal to noise ratio for such 
a received signal. 

These and other advantages of the present invention will become apparent 
upon reading the following detailed descriptions and studying the various 
10 figures of the drawings. 

Brief Description of the Drawings 

Figure 1 depicts a prior art receiver of a Multi-Carrier Modulated (MCM) signal 
sensed by antenna element 100 and fed 102 to module 110; 

Figure 2A depicts Frequency Division Multiple Access (FDMA) protocols, 
15 including GSM, using frequency bins separated by guard bands as the sub- 
carriers; 

Figure 2B depicts Orthogonal Frequency Division Multiplexing (OFDM) using 
densely spaced sub-carriers and overlapping spectra; 

Figure 2C depicts an OFDM technique used in conjunction with Direct 
20 Sequence Spread Spectrum techniques, to create CDMA-OFDM protocols; 
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Figure 2D depicts an example emitted power spectrum requirement for 
transmitted OFDM signals as found in Figure 120 of the IEEE 802.11a 
protocol specification; 

Figure 3 depicts a system 250 for processing a sample list 260 of at least two 
5 digitized samples based upon at least one electromagnetic receptor 100 using 
at least a truncated pseudo-inverse B1 410; 

Figure 4A depicts a detail flowchart of program system 1000 of Figure 3 
presenting program steps residing in memory 290 accessibly coupled 282 to 
processing engine 280; 

10 Figure 4B depicts an alternative flowchart of program system 1000 of Figure 3 
presenting program steps residing in memory 290 accessibly coupled 282 to 
processing engine 280; 

Figure 5A depicts a detail flowchart of operation 1012 of Figure 4A further 
receiving the sample list based upon the electromagnetic receptor; 

15 Figure 5B depicts a detail flowchart of operation 1112 of Figure 5A further 
receiving the first sample list based upon the first electromagnetic receptor; 

Figure 5C depicts a detail flowchart of operation 1122 of Figure 5A further 
receiving the second sample list based upon the second electromagnetic 
receptor; 

20 Figure 6A depicts a detail flowchart of operation 1022 of Figure 4A further 
processing the received sample list; 
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Figure 6B depicts a detail flowchart of operation 1022 of Figure 4A further 
processing the received sample list; 

Figure 6C depicts a matrix view of truncated pseudo-inverse B1 410 
comprised of third truncated pseudo-inverse B3 412 and a fourth truncated 
pseudo-inverse B4 41 4, with B3 412 vertically arranged with respect to B4 
414 and the effect of truncation upon these pseudo-inverses being the 
removal of columns; 

Figure 6D depicts a matrix view of truncated pseudo-inverse B1 410 
comprised of third truncated pseudo-inverse B3 412 and a fourth truncated 
pseudo-inverse B4 414, with B3 412 horizontally arranged with respect to B4 
414 and the effect of truncation upon these pseudo-inverses being the 
removal of rows; 

Figure 6E depicts a matrix view of truncated pseudo-inverse B1 410 
comprised of third truncated pseudo-inverse B3 412 and a fourth truncated 
pseudo-inverse B4 414, with B3 412 vertically and horizontally arranged with 
respect to B4 414 and the effect of truncation upon these pseudo-inverses 
being the removal of rows; 

Figure 7A depicts a detail flowchart of operation 1022 of Figure 4A further 
processing the received sample list; 

Figure 7B depicts a detail flowchart of program system 1000 of Figure 3 
further determining availability of a first communications service based upon 
the third truncated pseudo-inverse B3 and of a second communications 
service based upon the fourth truncated pseudo-inverse B4; 



29 



Figure 8 depicts a detail flowchart of program system 1000 of Figure 3 further 
determining availability of a first communications service based upon the third 
truncated pseudo-inverse B3 and of a second communications service based 
upon the fourth truncated pseudo-inverse B4 and using the method 
processing the sample list; 

Figure 9A depicts a detail flowchart of operation 1272 of Figure 7A further 
determining the first communications service availability; 

Figure 9B depicts a detail flowchart of operation 1282 of Figure 7A further 
determining the second communications service availability; 

Figure 10A depicts a detail flowchart of operation 1000 of Figures 7B and/or 8 
further performing the method of determining availability of a first and second 
communications service based upon the truncated pseudo-inverses; 

Figure 10B depicts a detail flowchart of program system 1000 of Figure 3 
determining of a Clear Channel Access for a CSMA protocol such as IEEE 
802.11a; and 

Figure 11 depicts an alternative embodiment of the invention from that of 
Figure 3 based upon at least one of a means for receiving the sample list 510 
and a means for using truncated pseudo-inverses 520; 

Figure 12 depicts simulation results comparing FFT to truncated DFT with a 
one vector constrained subspace processing sample lists for a 64-QAM 
modulation scheme such as IEEE 802.11a employs across a channel 
exhibiting AWGN with one known interferes 



Figure 13 depicts simulation results comparing FFT to truncated DFT with an 
eleven vector constrained subspace processing sample lists for a 64-QAM 
modulation scheme such as IEEE 802.11a employs across a channel 
exhibiting AWGN with one known interferes 

Detailed Description of the Invention 

Figure 3 depicts a system 250 for processing a sample list 260 of at least two 
digitized samples based upon at least one electromagnetic receptor 100 using 
at least one truncated pseudo-inverse B1 410. Multiple truncated pseudo- 
inverses, such as B2 420, may be preferred in certain embodiments of the 
invention. 

At least one processing engine 280 is receptively coupled 272-270-262-260- 
192-190-162-160-134-130-122-120-112-110-102 to the electromagnetic 
receptors 100 to provide the digitized samples 192. Processing engine 280 
may be preferably controlled by a program system 1000 comprising program 
steps residing in memory 290 accessibly coupled 282 to the processing 
engine 280. 

Program system 1000 preferably implements the inventive methods of 
operation. Alternatively, processing engine 280 may be hardwired to perform 
at least some of the steps of the methods described herein as implemented 
by program system 1000. 

The method of operation is thus not reliant upon a device resembling a 
computer. The inventive operations discussed herein may be embodied by a 
variety of means besides a computer. By way of example, systems employing 
one or a combination of at least one of program counter driven instruction 
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processor, finite state machines and pipelined dedicated processor engines, 
may be preferred for certain embodiments of the invention. The discussion of 
Figure 1 1 will further point out some examples of such embodiments of the 
invention. 

5 Processing engine 280 is also receptively coupled 272-270-262-260-212-210- 
172-170-136-130-122-120-112-110-102 to the electromagnetic receptors 100 
to provide the digitized samples 212. 

Figure 4A depicts a detail flowchart of program system 1000 of Figure 3 
presenting program steps residing in memory 290 accessibly coupled 282 to 
io processing engine 280. 

Arrow 1010 directs the flow of execution from starting operation 1000 to 
operation 1012. Operation 1012 performs receiving the sample list based 
upon the electromagnetic receptor to create a received sample list containing 
at least two received samples. Arrow 1014 directs execution from operation 
15 1012 to operation 1016. Operation 1016 terminates the operations of this 
flowchart. 

Arrow 1020 directs the flow of execution from starting operation 1000 to 
operation 1022. Operation 1022 performs processing the received sample list 
using the truncated pseudo-inverse B1 on at least some of the received 
20 samples to create a received signal parameter list containing at least one 
received signal parameter. Arrow 1024 directs execution from operation 1022 
to operation 1016. Operation 1016 terminates the operations of this flowchart. 

Note that the sample list is based upon a transported version of a signal 
progression generated using at least part of linear transform A1. The 
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transported version of a signal progression will usually be modified by noise, 
which can be from any of a number of mechanisms. Noise mechanisms may 
include, but are not limited to, background noise of the physical transport 
layer(s), thermal noise found in the reception mechanism typified by the upper 
5 half of Figure 3, as well as other potentially digital sources of noise including 
rounding errors in preceding calculations and digitization errors in the A/D 
converters. 

Certain embodiments of the invention view the transport mechanism itself as 
including a linear transform A1 acting upon the baseband time progression in 
10 terms of its reception at one or more receptors and/or one or more antenna 
elements of one or more receptors. 

Examples of such phenomena include, but are not limited to, radio signal 
propagation effects in air, often over time. These transforms can be derived 
from the propagation antenna lobe plots, which are usually done in polar 

15 coordinates. Rectangular grids in a polar coordinate system can be used to 
derive an angular versus temporal (propagation distance) based linear 
transform reasonably modeling the propagation effects. Such linear 
transforms are almost always singular. Often these linear transforms are 
between vector spaces of differing dimension, showing far more gradations of 

20 time/distance than angular gradations. Nonetheless, these transforms can be 
derived from the empirical lobe plots of the specific radio network. As such, 
their pseudo-inverses can be calculated and used to improve location 
resolution and/or improve the signal to noise ratio of the received signal. 

In some circumstances, receivers can be found to have significant non-linear 
25 effects. Common causes of such effects include, but are not limited to, the 
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non-linear effects of power transistors, amplifiers, among other effects. In 
certain circumstances, such non-linear effects can be adequately modeled 
using filter banks, which frequently employ a Discrete Wavelet Transform 
decomposition of the non-linear effect. When such a filter model can be 
5 accurately derived, a pseudo-inverse of the DWT matrix will yield the minimal 
least square estimate on the time varying stimulus driving the non-linear 
circuit element or system. 

Certain embodiments of the invention may include two electromagnetic 
receptors receptively coupled to processing engine 250. 

10 Figure 4B depicts an alternative flowchart of program system 1000 of Figure 3 
presenting program steps residing in memory 290 accessibly coupled 282 to 
processing engine 280. 

Arrow 1030 directs the flow of execution from starting operation 1000 to 
operation 1032. Operation 1032 performs receiving the sample list based 

is upon the electromagnetic receptor to create a received sample list containing 
at least two received samples. Arrow 1034 directs execution from operation 
1032 to operation 1036. Operation 1036 performs processing the received 
sample list using the truncated pseudo-inverse B1 on at least some of the 
received samples to create a received signal parameter list containing at least 

20 one received signal parameter. Arrow 1038 directs execution from operation 
1036 to operation 1040. Operation 1040 terminates the operations of this 
flowchart. 

One of skill in the art will recognize that Figures 4A and 4B depict equivalent 
operations, either of which may be preferred in different embodiments of the 
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invention. Figure 4A presents an execution mechanism commonly found in a 
reai-time, often event driven, operating environment. Figure 4B presents the 
operations as sequentially following one another, which is an approach often 
favored by some applications programming environments, emphasizing a flow 
5 of operations. In both Figures 4A and 4B, some or all of the depicted 
operations may be performed concurrently. 

Certain embodiments of the invention may preferably interact with different 
electromagnetic receptors. 

The electromagnetic receptor may include at least one antenna element, in 
10 which case the sample list based upon the electromagnetic receptor may be 
further derived from an electromagnetic field proximate with the antenna 
element. The electromagnetic receptor may include at least two antenna 
elements, in which case the sample list based upon the electromagnetic 
receptor may be further derived from the electromagnetic fields proximate with 
15 the antenna elements. 

The electromagnetic receptor may include at least one semiconductor 
receptor element, in which case the sample list based upon the 
electromagnetic receptor may be further derived from an electromagnetic field 
based upon the bulk transport properties of the semiconductor receptor 
20 element. The electromagnetic receptor may further include at least two 
semiconductor receptor elements, in which case the sample list based upon 
the electromagnetic receptor may be further derived from an electromagnetic 
field based upon the bulk transport properties of the semiconductor receptor 
elements. 



Examples of antenna elements include, but are not limited to, wire antennas, 
dipoles, quadrapoles, antenna arrays, horn antennas, and patch antennas, by 
way of example. Examples of semiconductor receptor elements include, but 
are not limited to, semiconductor lasers, masers, and Light Emitting Diodes 
(LEDs), by way of example. Semiconductor receptor elements may contain 
either crystalline materials and/or amorphous materials. Semiconductor 
receptor elements may contain either inorganic and/or organic compounds. 

The electromagnetic receptor may be comprised of a first electromagnetic 
receptor and a second electromagnetic receptor. Such embodiments include 
receivers coupled to multiple antenna sites, as well as receivers coupled to 
multiple semiconductor receptors. 

Figure 5A depicts a detail flowchart of operation 1012 of Figure 4A further 
receiving the sample list based upon the electromagnetic receptor, 

Arrow 1110 directs the flow of execution from starting operation 1012 to 
operation 1112. Operation 1112 performs receiving a first sample list based 
upon the first electromagnetic receptor to create the received sample list 
containing at least two received samples. Arrow 1114 directs execution from 
operation 1 1 1 2 to operation 1116. Operation 1116 terminates the operations 
of this flowchart. 

Arrow 1120 directs the flow of execution from starting operation 1012 to 
operation 1122. Operation 1122 performs receiving a second sample list 
based upon the second electromagnetic receptor to create the received 
sample list containing at least two received samples. Arrow 1124 directs 



execution from operation 1122 to operation 1116. Operation 1116 terminates 
the operations of this flowchart. 

Note that receiving the sample list based upon the electromagnetic receptor 
may include at least one of the two performed operations of Figure 5A. 

Figure 5B depicts a detail flowchart of operation 1112 of Figure 5A further 
receiving the first sample list based upon the first electromagnetic receptor. 

Arrow 1150 directs the flow of execution from starting operation 1112 to 
operation 1152. Operation 1152 performs receiving a first sample list based 
upon the first electromagnetic receptor to create a first received sample list 
containing at least two received first samples. Arrow 1154 directs execution 
from operation 1152 to operation 1156. Operation 1156 terminates the 
operations of this flowchart. 

Figure 5C depicts a detail flowchart of operation 1122 of Figure 5A further 
receiving the second sample list based upon the second electromagnetic 
receptor. 

Arrow 1170 directs the flow of execution from starting operation 1122 to 
operation 1172. Operation 1172 performs receiving a second sample list 
based upon the second electromagnetic receptor to create a second received 
sample list containing at least two received second samples. Arrow 1174 
directs execution from operation 1172 to operation 1176. Operation 1176 
terminates the operations of this flowchart. 

Figure 6A depicts a detail flowchart of operation 1022 of Figure 4A further 
processing the received sample list. 
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Arrow 1190 directs the flow of execution from starting operation 1022 to 
operation 1192. Operation 1192 performs processing the first received 
sample list by using the truncated pseudo-inverse B1 on at least some of the 
first received samples to create a first received signal parameter list 
5 containing at least one first received signal parameter. Arrow 1194 directs 
execution from operation 1192 to operation 1196. Operation 1196 terminates 
the operations of this flowchart. 

Arrow 1200 directs the flow of execution from starting operation 1022 to 
operation 1202. Operation 1202 performs processing the second received 
10 sample list by using the truncated pseudo-inverse B1 on at least some of the 
second received samples to create a second received signal parameter list 
containing at least one second received signal parameter. Arrow 1204 directs 
execution from operation 1202 to operation 1196. Operation 1196 terminates 
the operations of this flowchart. 

is Note that it may be preferable to include just one of the two performed 
operations of Figure 6A. 

Note that it may be preferable to use more than one pseudo-inverse on the 
second received samples as shown in Figure 3. 

Figure 6B depicts a detail flowchart of operation 1022 of Figure 4A further 
20 processing the received sample list. 

Arrow 1230 directs the flow of execution from starting operation 1022 to 
operation 1232. Operation 1232 performs processing the second received 
sample list by using a second truncated pseudo-inverse B2 420 of Figure 3 on 
at least some of the second received samples to create a second received 
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signal parameter list containing at least one second received signal 
parameter. Arrow 1234 directs execution from operation 1232 to operation 
1236. Operation 1236 terminates the operations of this flowchart. 

Truncated pseudo-inverse B1 410 of Figure 3 may preferably contain at least 
5 a third truncated pseudo-inverse B3 and a fourth truncated pseudo-inverse 
B4. 

Figure 6C depicts a matrix view of truncated pseudo-inverse B1 410 
comprised of third truncated pseudo-inverse B3 412 and a fourth truncated 
pseudo-inverse B4 414, with B3 412 vertically arranged with respect to B4 
io 414 and the effect of truncation upon these pseudo-inverses being the 
removal of columns. 

Figure 6D depicts a matrix view of truncated pseudo-inverse B1 410 
comprised of third truncated pseudo-inverse B3 412 and a fourth truncated 
pseudo-inverse B4 414, with B3 412 horizontally arranged with respect to B4 
15 414 and the effect of truncation upon these pseudo-inverses being the 
removal of rows. 

Figure 6E depicts a matrix view of truncated pseudo-inverse B1 410 
comprised of third truncated pseudo-inverse B3 412 and a fourth truncated 
pseudo-inverse B4 414, with B3 412 vertically and horizontally arranged with 
20 respect to B4 414 and the effect of truncation upon these pseudo-inverses 
being the removal of rows. 

Figures 6C through 6E depict some embodiments of a truncated pseudo- 
inverse composed of more than one truncated pseudo-inverse. Such 
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compositions may include pseudo-inverses of different types, B3 may be a 
Drazin pseudo-inverse and B4 may be a Moore-Penrose inverse, for example. 

Further note that the truncation process may remove rows or columns which 
are not adjacent to each other. 

5 Figure 7 A depicts a detail flowchart of operation 1022 of Figure 4A further 
processing the received sample list. 

Arrow 1250 directs the flow of execution from starting operation 1022 to 
operation 1252. Operation 1252 performs processing the second received 
sample list by using the third truncated pseudo-inverse B3 on at least some of 
10 the received samples to create a third received signal parameter list 
containing at least one third received signal parameter. Arrow 1254 directs 
execution from operation 1252 to operation 1256. Operation 1256 terminates 
the operations of this flowchart. 

Arrow 1260 directs the flow of execution from starting operation 1022 to 
15 operation 1262. Operation 1262 performs processing the second received 
sample list by using the fourth truncated pseudo-inverse B3 on at least some 
of the received samples to create a fourth received signal parameter list 
containing at least one fourth received signal parameter. Arrow 1264 directs 
execution from operation 1262 to operation 1256. Operation 1256 terminates 
20 the operations of this flowchart. 

Certain embodiments of the invention include a method of determining 
availability of a first communications service based upon the third truncated 
pseudo-inverse B3 and of a second communications service based upon the 
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fourth truncated pseudo-inverse B4 using the method processing the sample 
list. 

Figure 7B depicts a detail flowchart of program system 1000 of Figure 3 
further determining availability of a first communications service based upon 
5 the third truncated pseudo-inverse B3 and of a second communications 
service based upon the fourth truncated pseudo-inverse B4. 

Arrow 1270 directs the flow of execution from starting operation 1000 to 
operation 1272. Operation 1272 performs determining the first 
communications service availability based upon the third received signal 
10 parameter list to create a first communication service determination. Arrow 
1274 directs execution from operation 1272 to operation 1276. Operation 
1276 terminates the operations of this flowchart. 

Arrow 1280 directs the flow of execution from starting operation 1000 to 
operation 1282. Operation 1282 performs determining the second 
15 communications service availability based upon the fourth received signal 
parameter list to create a second communication service determination. 
Arrow 1284 directs execution from operation 1282 to operation 1276. 
Operation 1276 terminates the operations of this flowchart. 

This method of determining availability of a first communications service and 
20 of a second communications service based upon the truncated pseudo- 
inverses and using the method processing the sample list may also be seen 
as a standalone application as follows. Such an application may preferably 
run upon wake-up in a transceiver. 
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Figure 8 depicts a detail flowchart of program system 1000 of Figure 3 further 
determining availability of a first communications service based upon the third 
truncated pseudo-inverse B3 and of a second communications service based 
upon the fourth truncated pseudo-inverse B4 and using the method 
5 processing the sample list. 

Arrow 1310 directs the flow of execution from starting operation 1000 to 
operation 1312. Operation 1312 performs receiving the sample list based 
upon the electromagnetic receptor to create a received sample list containing 
at least two received samples. Arrow 1314 directs execution from operation 
10 1312 to operation 1316. Operation 1316 terminates the operations of this 
flowchart. 

Arrow 1320 directs the flow of execution from starting operation 1000 to 
operation 1322. Operation 1322 performs processing the received sample list 
by using the third truncated pseudo-inverse B3 on at least some of the 
15 received samples to create a third received signal parameter list containing at 
least one third received signal parameter. Arrow 1324 directs execution from 
operation 1322 to operation 1316. Operation 1316 terminates the operations 
of this flowchart. 

Arrow 1330 directs the flow of execution from starting operation 1000 to 
20 operation 1332. Operation 1332 performs processing the received sample list 
by using the fourth truncated pseudo-inverse B4 on at least some of the 
received samples to create a fourth received signal parameter list containing 
at least one fourth received signal parameter. Arrow 1334 directs execution 
from operation 1332 to operation 1316. Operation 1316 terminates the 
25 operations of this flowchart. 
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Arrow 1340 directs the flow of execution from starting operation 1000 to 
operation 1342. Operation 1342 performs determining the first 
communications service availability based upon the third received signal 
parameter list to create a first communication service determination. Arrow 
1344 directs execution from operation 1342 to operation 1316. Operation 
1316 terminates the operations of this flowchart. 

Arrow 1350 directs the flow of execution from starting operation 1000 to 
operation 1352. Operation 1352 performs determining the second 
communications service availability based upon the fourth received signal 
parameter list to create a second communication service determination. 
Arrow 1354 directs execution from operation 1352 to operation 1316. 
Operation 1316 terminates the operations of this flowchart. 

Figure 9A depicts a detail flowchart of operation 1272 of Figure 7 A further 
determining the first communications service availability. 

Arrow 1390 directs the flow of execution from starting operation 1272 to 
operation 1392. Operation 1392 performs detecting system communication 
based upon the third received signal parameter list to create a first system 
channel detection. Arrow 1394 directs execution from operation 1392 to 
operation 1396. Operation 1396 performs generating the first 
communications service determination based upon the first system channel 
detection. Arrow 1398 directs execution from operation 1396 to operation 
2400. Operation 2400 terminates the operations of this flowchart. 

Figure 9A depicts the determination of the availability of communications 
service for a communication protocol relying upon at least one system 
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channel being able to be detected and decoded. Examples of such MCM 
communications protocols include, but are not limited to, AMPs, GSM, 1S-95, 
Edge, and W-CDMA. 

Figure 9B depicts a detail flowchart of operation 1282 of Figure 7 A further 
5 determining the second communications service availability. 

Arrow 1410 directs the flow of execution from starting operation 1282 to 
operation 1412. Operation 1412 performs estimating a first energy term 
based upon the fourth received signal parameter list and a second energy 
term based upon the fourth received signal parameter list. Arrow 1414 directs 
10 execution from operation 1412 to operation 1416. Operation 1416 performs 
generating the second communications service determination based upon the 
first energy estimate exceeding the second energy estimate multiplied by a 
threshold value. Arrow 1418 directs execution from operation 1416 to 
operation 1420. Operation 1420 terminates the operations of this flowchart. 

15 Certain embodiments of the invention determine communications service 
availability by calculating a Clear Channel Assessment, CCA_flag for the 
communications protocol by estimating at least two energy terms based upon 
the value list. And generating the second communications service 
determination based upon whether the first energy term exceeds the second 

20 energy term multiplied by a threshold. 

Such embodiments do not have to wait for synchronization with a training 
sequence or preamble. They rely instead upon the physical characteristics of 
the encoded channel of the communications protocol, which must expend 
energy above the noise floor for the signal to be received. 
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Such embodiments are applicable to the OFDM communications protocols in 
general, and to communication protocols compatible with the IEEE 802.11 
specification in particular. 

Some embodiments of the invention use estimates of the peak power versus 
5 the average power for the two energy terms, while other embodiments 
estimate the channel signal energy and the channel noise energy as the two 
energy terms. 

Note that in certain embodiments of the invention, it may be preferred that 
both operation 1272 and 1282 determining distinct communication service 
io capabilities employ similar mechanisms, that is, both may employ only one of 
the mechanisms of Figures 9A and 9B. 

One of skill in the art will see that while the performed operations of Figures 
9A and 9B are shown in an essentially sequential flow of control, they may 
equivalently be implemented in a concurrent real-time operating paradigm. 
15 The choice of portrayal in Figures 9A and 9B was made strictly to clarify the 
discourse and is not meant to limit the scope of the claims. 

In certain embodiments of the invention the second communications service 
determination may preferably include a Clear Channel Access determination 
in a fashion applicable to at lease CSMA communications protocols. 

20 The linear transform A may preferably include an FFT. Linear transform A 
may preferably be the FFT of 64 points as specified in IEEE 802.11a. 
Truncated pseudo-inverse B1 may preferably provide at least an 
approximation of the 52 active frequency bins of the IEEE 802.11a physical 
layer during data transmission. 



Figure 10A depicts a detail flowchart of operation 1000 of Figures 7B and/or 8 
further performing the method of determining availability of a first and second 
communications service based upon the truncated pseudo-inverses. 

Arrow 1450 directs the flow of execution from starting operation 1000 to 
5 operation 1452. Operation 1452 performs processing the received sample list 
by using a fifth truncated pseudo-inverse B5 on at least some of the received 
samples to create a second received signal parameter list containing at least 
one second received signal parameter. Arrow 1454 directs execution from 
operation 1452 to operation 1456. Operation 1456 terminates the operations 
io of this flowchart. 

The fifth truncated pseudo-inverse B5 may preferably provide at least an 
approximation of the 12 active frequency bins of the IEEE 802.11a physical 
layer during header transmission. 

Note that the operations described through Figure 9B further describe what 
15 may preferably be a standalone application in certain embodiments of the 
invention supporting determination of a Clear Channel Access for a CSMA 
protocol such as IEEE 802.1 1a. 

Figure 10B depicts a detail flowchart of program system 1000 of Figure 3 
determining of a Clear Channel Access for a CSMA protocol such as IEEE 
20 802.11a. 

Arrow 1470 directs the flow of execution from starting operation 1000 to 
operation 1472. Operation 1472 performs receiving the sample list based 
upon the electromagnetic receptor to create a received sample list containing 
at least two received samples. Arrow 1474 directs execution from operation 
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1472 to operation 1476. Operation 1476 performs processing the received 
sample list by using the truncated pseudo-inverse B1 to at least some of the 
received samples to create a received signal parameter list containing at least 
one received signal parameter. Arrow 1478 directs execution from operation 

5 1476 to operation 1480. Operation 1480 performs estimating a first energy 
term based upon the received signal parameter list and a second energy term 
based upon the received signal parameter list. Arrow 1482 directs execution 
from operation 1480 to operation 1484. Operation 1484 performs generating 
the communications service determination based upon the first energy 

10 estimate exceeding the second energy estimate multiplied by a threshold 
value. Arrow 1486 directs execution from operation 1484 to operation 1488. 
Operation 1488 terminates the operations of this flowchart. 

One of skill in the art will recognize that the sample list may be further based 
upon a transported version of a baseband signal progression generated using 
15 at least part of linear transform A1 . 

The baseband signal progression may be further generated as a signal 
progression using at least part of a second linear transform A2. 

The second linear transform A2 may be approximately similar to a Hadamard 
transform. Second linear transform A2 may be further approximately similar to 
20 a Walsh-Hadamard transform. 

The truncating pseudo-inverse B1 may approximate a pseudo-inverse of at 
least part of linear transform A1 applied to at least part of a second linear 
transform A2. 
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Certain embodiments of the invention may include at least part of the second 
linear transform A2 providing a scattering transform applied to a time 
progression generated by using at least part of a spreading linear transform 

A3. 

5 The truncating pseudo-inverse B1 may approximate a pseudo-inverse of at 
least part of linear transform A1 applied to at least part of the second linear 
transform A2 applied to at least part of the spreading linear transform A3. 

Figure 11 depicts an alternative embodiment of the invention from that of 
Figure 3 based upon at least one of a means for receiving the sample list 510 
10 and a means for using truncated pseudo-inverses 520. 

Processor 500 may embody at least some of the steps of the inventive 
methods as separate means for performing those steps. 

At least one means 510 is receptively coupled 274-270-262-260-192-190- 
162-160-134-130-122-120-112-110-102 to the electromagnetic receptors 100 
15 to provide the digitized samples 192. 

The inventive operations discussed herein may be embodied by a variety of 
means besides a computer. By way of example, systems employing one or a 
combination of at least one of the following, program counter driven 
instruction processing, finite state machines and pipelined dedicated 
20 processor engines may be preferred for certain embodiments of the invention. 

Means 510 may also receptively coupled 274-270-262-260-212-210-172-170- 
136-130-122-120-112-110-102 to the electromagnetic receptors 100 to 
provide the digitized samples 212. 
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Means 510 performs at least the operations of step 1012 of Figure 4 A. This 
includes, but is not limited to, receiving the sample list based upon the 
electromagnetic receptor to create a received sample list 400 containing at 
least two received samples. 

5 Means 510 may employ one or a combination of mechanisms including, but 
not limited to, at least program counter driven instruction processing, finite 
state machines and pipelined dedicated processor engines for certain 
embodiments of the invention. 

Means 510 provides a mechanism by which the digitized samples become 
10 512 the received sample list 400, which can be readably accessed 522 by 
means 520. Note that it may be preferable in certain embodiments of the 
invention that means 520 be able to assert signals 524 for received sample 
list 400. Such asserted signals 524 may include, but are not limited to, 
addressing and control signals regulating which received samples or 
15 components of received samples are to be read, as well as potentially the 
ability to write data to the received sample list 400. Such embodiments may 
preferably support in-place calculations similar to in-place FFT calculations. 

Means 520 performs at least the operations of step 1022 of Figure 4A. This 
includes, but is not limited to, processing the received sample list by using the 
20 truncated pseudo-inverse B1 410 on at least some of the received samples in 
received sample list 400 to create a received signal parameter list containing 
at least one received signal parameter. 

Means 520 may employ one or a combination of mechanisms including, but 
not limited to, at least program counter driven instruction processing, finite 
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state machines and pipelined dedicated processor engines for certain 
embodiments of the invention. 

Consider the following simulation experiment. Simulate a simple TX-channel- 
RX simulation chain (based on IEEE 802.11a specification) featuring two 

5 receivers: the first receiver is FFT-based receiver while the second one is a 
constrained DFT-based received. The transmitter is an IFFT-based 
transmitter compliant with the IEEE 802.11a specification. The channel is 
AWGN (subject to Average White Gaussian Noise) and a known interferer U 
is added to the transmit signal. The TX data- and pilot-bearing subcarriers (48 

io + 4 = 52 subcarriers referred to as "a") are saved and used in conjunction with 
the RX ones (referred to as "a") to construct a figure of merit which is the 
following normalized dot product: 

FOM = ( a - a ) * { a - a )' / ( a * a' ) 

where superscript 1 designates the hermitian transpose operator. The 
is denominator is a normalizing quantity (that has the dimension of energy like 
the numerator). 

The FOM is an indicator of the degree of proximity of frequency-domain 
vectors a and a. Ideally, in the absence of any noise or interferer or 
imperfections of any sort, a = a and FOM - 0. The smaller the FOM the 
20 better. 

The flow of operations is as follows: 

1. Frequency-domain signal "a" is transformed by the IFFT to produce time- 
domain signal x (TX) 
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2. Noise as well as single interferer U are added to x to produce time-domain 
signal y = x + U + scaled noise 

3. Time-domain signal y is transformed by an FFT to produce frequency- 
domain signal a (RX) 

5 4. FOM is computed 

The constrained DFT matrix was generated as follows using Matlab: if M is 
the 64 x 52 DFT matrix obtained by truncating the full 64 x 64 matrix 
associated with I FFT and U is the 64-element column vector associated with 
the undesired interferer to be removed then the constrained DFT matrix to be 
used in the receiver is given by: 

N = pinv([M U] ) ; % compute pseudo-inverse (U may also 
designate a collection of vectors rather than a 
single one) 

N = N(l:52,:); % truncate 53 (or whatever) x 64 

pseudo-inverse to obtain a 52 x 64 matrix 
a = N*y; % perform demodulation 

The simulation results clearly show the effect of using the constrained DFT 
over the FFT. By removing the expected interferer, the use of N reduces the 
gap between transmit and receive modulated signals (FOM converges 
towards zero or a noise-dependant threshold) thereby dramatically reducing 
the Burst Error Rate (BER) and/or Peak Error Rate(PER) loss due to U. 

The residual FOM is due to the presence of AWGN, whose contribution 
cannot be eliminated by either the constrained DFT or by the FFT. The 
constrained DFT cancels only those signals that are a linear combination of 
25 the vectors spanning the subspace U. 
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Figure 12 depicts simulation results comparing FFT to truncated DFT with a 
one vector constrained subspace processing sample lists for a 64-QAM 
modulation scheme such as IEEE 802.11a employs across a channel 
exhibiting AWGN with one known interferes 

5 Figure 13 depicts simulation results comparing FFT to truncated DFT with an 
eleven vector constrained subspace processing sample lists for a 64-QAM 
modulation scheme such as IEEE 802.11a employs across a channel 
exhibiting AWGN with one known interferes 

In both Figure 12 and 13 the vertical axis represents the normalized Root 
10 Mean Square (RMS) error between a and a. This is the figure of merit in 
comparing implementations of differing constraint subspace dimension. 

In both Figure 12 and 13 the horizontal axis represents the Signal to Noise 
Ratio (SNR) in decibels (dB). 

The simulation results indicate that the immunity to AWGN deteriorates as the 
15 dimension of the constraining subspace increases. There is a tradeoff 
between immunity to AWGN and the number of interferers (or to be more 
precise the dimension of the constraining subspace U) that can be cancelled. 
The simulation results show a penalty of a low SNR the larger the dimension 
of the constraint subspace. 

20 Please note that it is also possible to remove any interferer from the received 
signal y by simply projecting y on the subspace orthogonal to U rather than 
using N. In this case, construct the projection operator as follows: 

P = I - U * inv(U l * U) * U ' 



This operator can be seen to be idempotent. Then, z = P * y is interferer-free, 
but can not be subsequently Fast Fourier Transformed since z is no longer 
colinear to x, because a and a are very likely to be far apart. 

The preceding embodiments have been provided by way of example and are 
5 not meant to constrain the scope of the following claims. 
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